On an efficient numerical method for modeling sea ice dynamics
نویسنده
چکیده
A computationally efficient numerical method for the solution of nonlinear sea ice dynamics models employing viscous-plastic rheologies is presented. The method is based on a semi-implicit decoupling of the x and y ice momentum equations into a form having better convergence properties than the coupled equations. While this decoupled form also speeds up solutions employing point relaxation methods, a line successive overrelaxation technique combined with a tridiagonal matrix solver procedure was found to converge particularly rapidly. The procedure is also applicable to the ice dynamics equations in orthogonal curvilinear coordinates which are given in explicit form for the special case of spherical coordinates.
منابع مشابه
Modeling Arctic sea ice with an efficient plastic solution
A computationally efficient numerical method is developed for solving sea ice momentum equations that employ a nonlinear viscous-plastic rheology. The method is based on an alternating direction implicit (ADI) technique that involves a direct solution of the momentum equations. This method is therefore more computationally efficient than those employing an iterative procedure in solving the equ...
متن کاملInelastic Continuum Modeling of Carbon Nanotube,s Behavior Using Finite Element Method
This paper describes a continuum model for analyzing the inelastic behavior of a single walled carbon nanotube (SWCNT) in different loading conditions. Because of limitations in using molecular dynamics (and other atomic methods) to model the failure load of the SWCNT, continuum mechanics methods are considered in this paper. Based on some experimental and theoretical results, an elasto-plastic...
متن کاملUsing the material-point method to model sea ice dynamics
[1] The material-point method (MPM) is a numerical method for continuum mechanics that combines the best aspects of Lagrangian and Eulerian discretizations. The material points provide a Lagrangian description of the ice that models convection naturally. Thus properties such as ice thickness and compactness are computed in a Lagrangian frame and do not suffer from errors associated with Euleria...
متن کاملAn Optimization Approach to Modeling Sea Ice Dynamics, Part 2: Finite Ice Strength Effects
The effects of a finite ice strength on a new model for sea ice dynamics, deriving the internal pressure field from a global optimization problem, rather than a local rheology, are examined. Building on the promising results from the one-dimensional Lagrangian model described previously, here we add one of the key properties of sea ice. In order to investigate the behavior of the model under ic...
متن کاملA numerical model for electromagnetic scattering from sea ice
A numerical model for scattering from sea ice based on the finite difference time domain (FDTD) technique is presented. The sea ice medium is modeled as consisting of randomly located spherical brine scatterers with a specified fractional volume, and the medium is modeled both with and without a randomly rough boundary to study the relative effects of volume and surface scattering. A Monte Carl...
متن کامل